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Abstract 

We investigate the stabilizing effects of the magnetic fields in the linearized magnetic Rayleigh- 
Taylor (RT) problem of a nonhomogeneous incompressible viscous magnetohydrodynamic fluid 
of zero resistivity in the presence of a uniform gravitational held in a three-dimensional bounded 
domain, in which the velocity of the fluid is non-slip on the boundary. By adapting a modi¬ 
fied variational method and careful deriving a priori estimates, we establish a criterion for the 
instability/stability of the linearized problem around a magnetic RT equilibrium state. In the 
criterion, we find a new phenomenon that a sufficiently strong horizontal magnetic field has the 
same stabilizing effect as that of the vertical magnetic field on growth of the magnetic RT in¬ 
stability. In addition, we further study the corresponding compressible case, i.e., the Parker (or 
magnetic buoyancy) problem, for which the strength of a horizontal magnetic field decreases with 
height, and also show the stabilizing effect of a sufficiently large magnetic field. 

Keywords: Viscous MHD equations; equilibrium state; magnetic Rayleigh-Taylor instability; 
incompressible/compressible fluids; Parker’s instability; instability/stability criteria. 
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1. Introduction 


The three-dimensional (3D) nonhomogeneous, incompressible and viscous magnetohydrody¬ 
namic (MHD) equations with zero resistivity (i.e. without magnetic diffusivity) in the presence of 
a uniform gravitational field in a domain D C reads as follows (see, for example, js], @, 29 

on the derivation of the equations): 


^ pt + div(pn) = 0, 

pvt + pv ■ Vv -I- V(p -I- AoV|Mp/2) = /iAn -1 - AqM ■ VM — pf/ea, 
Mt = M ■ Vn — V ■ VM, 
divn = divM = 0. 


Here the unknowns p = p{x,t), v = v{x,t), M := M{x,t) and p = p{x,t) denote the density, 
velocity, magnetic field and pressure of the incompressible fluid, respectively; p > 0 stands for 
the coefficient of shear viscosity, Aq for the permeability of vacuum, g > 0 for the gravitational 
constant, 63 = (0, 0,1)"*" for the vertical unit vector, and —pge^ for the gravitational force. In 
the system fll.ip the equation fll.lD i is the continuity equation, fll.lD o describes the balance law 
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of momentum, while fll.ip -i is called the induction equation. As for the constraint divM = 0, 
it can be seen just as a restriction on the initial value of M since (divM)^ = 0 due to ( 1 . 1 ) 3 . 
We remark that the resistivity is neglected in fll.ip q. and this arises in the physics regime with 
negligible electrical resistance. 

Now, we choose a Rayleigh-Taylor (RT) density prohle p := p^x^), which is independent of 
{xi,X 2 ) and satishes 


peC^(n), 


inf p > 0, 

xen 


P 


X3=X^ 


> 0 


for some x^ G {^3 | {xi,X 2 ,X 3 ) G fl}. 


( 1 , 2 ) 


where p' := dp/da; 3 , x^ is the third component of xq G fl. We refer to 2l|, Remark 1.1] for 


the construction of such p. We remark that the second condition in fll.2p prevents us from 
treating vacuum in the construction of unstable solutions, while the third one in fll.2p assures 
that there is at least a region in which the RT density has larger density with increasing height 
X 3 and may lead to the classical RT instability. Since we investigate the stabilizing effects of the 
horizontal and vertical magnetic helds, we consider the magnetic held prohle M = mei (i.e. the 
magnetic held orthogonal to the direction of the gravitational force) or me^, (i.e., the magnetic 
held parallel to the direction of the gravitational force), where Ci := (1, 0, 0)"*" and m is non-zero 
constant. Then (p, u, M) = (p, 0, M) dehnes a magnetic RT equilibrium state to fll.ip . where the 
equilibrium pressure prohle p is determined by 


Vp = -ppe3. 

Denoting the perturbation to the RT equilibrium state by 

Q = p — p, u = V — 0, N = M — M, q = P — P, 
then, {g,u,q) satishes the perturbed equations 

Pt + M • V(p-hp) = 0, 

(p + p)ut + (p + p)u ■ Vtt -|- V(p' -|- AqI+ ii^p/2) 

< = pAu +Xo{N + M)-ViN + M) - {g +p)ge 3 , (I. 3 ) 

Nt = {N + M) -Vu-u- V{N + M), 
divM = divA^ = 0. 

For the system fll.Sp . we impose the initial and boundary conditions: 

(p, u, N)\t=o = (po. Mo, No) in D, (1.4) 

u{x,t)\gfi = 0 for any t > 0. (1.5) 

If the perturbation {g,u,N) is very small, then the two-order small terms (i.e., the nonlinear 
terms) could be neglected in the perturbed equations, and we obtain the following linearized 
magnetic RT equations around the equilibrium state (p, 0,M): 

Qt + p'u3 = 0 , 

put + V{q + XomNi) = pAu -F XomdiN - gges, 

AT O 

JXIt = mOiU, 

divM = divA^ = 0, i = 1,3. 
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The linearized equations are convenient in mathematical analysis in order to have an insight 
into the physical and mathematical mechanisms of the magnetic RT instability. In this article, 
we investigate the instability/stability of the linearized magnetic RT problem fll.4|) - fll.6l) in a 
bounded domain. 

Next we briefly introduce the related background and research motivation of the magnetic 
RT instability. The RT instability is well-known as a gravity-driven instability in fluid dynamics 
when a heavy fluid is on top of a light one. This phenomenon was hrst studied by Rayleigh 1883 
and then Taylor, thus called Rayleigh-Taylor instability. The analogue of the RT instability 


arises when fluids are electrically conducting and a magnetic held is present, and the growth 
of the instability will be inhuenced by the magnetic held due to the generated electromagnetic 
induction and the Lorentz force. In the last decades, the stabilizing ehect of the magnetic held on 
RT instability has been analyzed by a number of authors. Kruskal and Schwarzchild in 1954 hrst 
showed that a horizontal magnetic held has no ehect on growth of the linear RT instability [^ . 
and such RT instability arising in magnetohydrodynamics is called the Kruskal-Schwarzschild 
instability. Then the stabilizing role of a vertical magnetic held was further investigated by 
Hide in 17|] where the ehect of hnite viscosity and resistivity was included and his analysis was 
encumbered with many parameters. These results were summarized in the monograph by 
Chandrasekhar. Similar results also hold for the problem fll.4p - fll.6l) dehned in a horizontally 
periodic domain (i.e., := (27rLT)^ x (—/,/), where 0 < / < -l-cxo and 27rLT stands for the 

ID-torus of length 2ttL). More precisely, one can construct an unstable solution to fll.4p - fll.6l) 
for any horizontal magnetic held M = mei and for the vertical magnetic held M = me^ with 
m < Me, where Mq denotes a critical number dehned by 


Me : = 


g f_ip'\'ip{x3)\Mx3 

\ - 1 - 

\ Aof_i \ij'{x3)\'^dX3 


> 0 , 


(1.7) 


see [22|, Theorem 1.1] for more details. Moreover, if I G (0,-|-cxd), one can show the stability of 
fll.4l) - fll.6l) for any m > Me, please refer to 37| for the proof of the linear stability. This means 
that a sufficiently strong vertical magnetic held can prevent the growth of the linear RT instability. 
These diherent ehects on instability between the horizonal and vertical magnetic helds motivate 
us to study the intrinsic mechanism from the mathematical point of view. By a simple analysis 
of the Lorentz force, we hnd that the non-slip boundary condition of the velocity can play an 
important role in the stabilizing ehect of the vertical magnetic held. On the other hand, we also 
notice that the velocity is horizontally periodic and non-slip at X 3 = ±/ for I G (0, -|-oo), i.e, the 
boundary conditions for the velocity are diherent in the horizontal and vertical directions. This 
shows that the diherent behaviors on instability of the horizonal and vertical magnetic helds may 
be somehow related to the boundary condition of the velocity. The same situation occurs for the 
inviscid magnetic RT problem [l^ (i.e. /i = 0 in fll.5l) - fll.6lB in a 2D periodic domain, where the 
velocity is periodic in both vertical and horizontal directions. A natural question arises whether 
the horizontal magnetic held has the same stabilizing ehect as the vertical one on growth of the 
magnetic RT instability, when the velocity is non-slip on the boundary of a bounded domain D. 

The hrst aim of this article is to give a positive answer, namely, we shall hnd that there is a 
critical number for a horizontal magnetic held M = me^, such that the linearized magnetic RT 
steady state is stable provided m > (see Theorem 12.ip . Moreover, we hnd that the horizontal 
magnetic held has the same stabilizing ehect as the vertical one on growth of the magnetic RT 
instability (see Remark [ 2 . 2 p . Thus, this article updates Kruskal and Schwarzschild’s results. To 
our best knowledge, this is the hrst article to study the role of the non-slip boundary condition 
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of velocity in the horizontal direction in the MHD instability under a horizontal magnetic field. 
We also mention that the magnetic RT instability for two-layer incompressible fluids separated 
by a free interface (stratihed fluids) in a horizontally periodic domain is investigated and some 
weakly nonlinear instability results are obtained, see |23| . 

The second aim is to extend Theorem 12.11 for the incompressible MHD problem fll.5p - fll.6p 
with M = mei to the corresponding compressible isentropic case, in which the horizontal mag¬ 
netic held is vertically stratihed. Before stating our result, we introduce the governing equations. 
The corresponding compressible isentropic model of fll.ip reads as follows. 

pt + div(pn) = 0 , 

pvt + pv-Vv + V{p + Xo\M\^/2)=pAv + poV(\.Y>JV + XQM-VM — pge^, 

Mt = M ■ Vv — V ■ VM — Mdivn, 

^ divM = 0 , 

where po ■.= p + i>, i> denotes the bulk viscosity and Su + 2p > 0. The pressure p is usually 
determined through the equations of state. In this article we focus our study on the isentropic 
how case and consider that 

P = p{p) = Ap^^ ( 1 - 9 ) 

where 7 > 1 denotes the adiabatic constant and H > 0 is a constant. 

Next we construct a magnetic RT equilibrium state for fll.Sp . Letting p satisfy fll.2l) . we dehne 
a horizontal magnetic held prohle Me := (me, 0 , 0 ) with 


rric = ±-1 / — {C — p{p) — gF{p)), (a function of X 3 ) 

V Ao 

where F{p) denotes a primitive function p and (7 is a positive constant satisfying 

C — p{p) — gF{p) >0 on D. 


( 1 . 10 ) 


It is easy to see that 

P'{p)p' = -Xorricm'^ - gp, (1.11) 

where p'{p) := Thus {p,0,Mc) constructed above is an equilibrium state to fll.Sp . i.e., 

V(p(p)-|-Ao|Mc|^/2) = AqMc • VMc — ppes and divMc = 0. (1.12) 

Obviously, by virtue of the relation fll.lip . rric is impossible to be a constant. We remark here 
that, by virtue of the relation fll.l2p . there does not exist a magnetic RT equilibrium state 
(p, 0, M) of fll.81) . such that M is a vertical magnetic held. 

Now denoting the perturbation to (p, 0, Me) by 

g = p — p, u = V — 0, N = M — Me, 


we get the perturbed equations: 

Pt-hdiv((p-Fp)M) = 0, 

{g + p)ut + {g + p)u ■ Vu + V{p{g + p) + Xo\N_+ /2) 

< = pAu +poVdivu +Xo{N + Me)-ViN + Me) - {g +p)ge 3 , (1T3) 

Nt = {N + Me) -Vu-u - V{N + Me) -{N + Me)divu, 
divA^ = 0. 
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We impose the following initial and boundary conditions for fll.ldp : 


{g,u,N)\t=o = {qo,Uo,No) inf], 
u{x,t)\gQ = 0 for any f > 0. 

The linearized equations of (ll.ldh around the equilibrium state (p, 0, M) read as 

Qt + diy{pu) = 0, 

put + V(p'(p)p + XorricNi) = pAu + poVdivM + AoiVsM' + XorricdiN - gge^, 
Nt = rricdiu — — Mcdivw, 

divA^ = 0, 


(1.14) 

(1.15) 


(1.16) 


where := dMc/dx^. The system fll.l6p with suitable initial and boundary conditions consti¬ 
tutes a compressible (viscous) magnetic RT problem. 

Now, we introduce the related research background on the compressible magnetic RT problem. 
By virtue of the conditions fll.2p . fll.9|) and fll.lip . there exists at least a region in which the RT 
density prohle has larger density with increasing height and the magnetic held Me causes a 
non-zero Lorentz force Xo{Mc ■ VM^c — V|McP/2), the direction of which is opposite to gravity, 
to support the heavier gas layered on top of lighter one. In particular, if the strength of the 
Lorentz force increases in fll.l2p . we see that the Lorentz force plays a role of buoyancy to drive 
the heavier gas to go up. This idea of “magnetic buoyancy” was introduced by Parker in 
connection with the formation of sunspots. When sunspots break out, most solar prominences 
will go up due to the magnetic buoyancy, and then slowly fall down toward the surface of the 
sun due to gravity of the sun, while some solar prominences can hoat in the solar corona for a 
long time, although the density of the former is heavier 1000-10000 times than the latter. Hence, 
if the equilibrium state fll.121) is slightly disturbed, then the instability caused by magnetic 
buoyancy and gravity may occur in compressible MHD problems. Therefore, such an instability 
is commonly referred as the magnetic buoyancy instability due to the new characteristic caused 
by the magnetic buoyancy, also called the Parker instability in the astronomical literature [35l |. 
Thus, the linearized compressible magnetic RT problem fll.14p - fll.16p can be called the linearized 
Parker problem, we refer the reader to [H, 24] and the references cited therein for more physical 
background on the Parker instability. The Parker instability is also suggested to be responsible for 
other observed astrophysical effects, for examples, Parker [35| demonstrated that the interstellar 
medium is unstable due to the magnetic buoyancy as an instability mechanism, and thought that 
the Parker instability is associated with interstellar cloud formation. It is worth pointing out 
that Fukui, et ah [sl observed giant molecular loops in a Galactic center, which offers a evidence 
for magnetic floating and supports thus Parker’s theory. 

The linear Parker instability have been widely investigated by physicists from the physical 
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and numerical simulation points of view, see [3, [if 
for example. In this paper, we study the dynamical instability and stability of the linearized 
Parker problem fll.14p - fll.16p from the mathematical point of view. More precisely, we give 
sufficiently conditions for the linear instability and stability of the problem fll.14p - fll.16p in the 
Hadamard sense in Sobolev spaces. Moreover, from the stability condition, we can easily see 
that a sufficiently strong Me has a remarkable stabilizing effect in the development of the Parker 
instability. Therefore, our results for the compressible can be regarded as a generalization of 
those for the previous incompressible case. 

The rest of this paper is organized as follows. In Section [2] we introduce the instability and 
stability of both linearized magnetic RT and Parker problems. Sections [3] and H] are devoted to 
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the proof of instability and stability of the two linearized problems. Finally, we give an additional 
resnlt on the critical nnmber in a horizontally periodic domain and prove the sharp growth rate 
of solntions to the two linearized problems in Section |5l 

2. Main results 

Before stating the main resnlts, we introdnce the notations nsed thronghont this paper. We 
always assnme that f] be a C^-smooth bonnded domain. For simplicity, we drop the domain f] 
in Sobolev spaces and the corresponding norms as well as in integrands over ff, for example, 

:= L\Q), Hi := := and 

We denote 

:= (0, cx)). Hi := {u E Hi \ divn = 0}, Jp) := J pw'^dx, 

A = {w E Hi I J{w) = 1} , Aa = An Hi, 
fj denotes the j-th component of the vector fnnction /q. 

The letter C denotes a generic positive constant which may depend on hi and other known 
physical qnantities snch as g, p, Xq, p, m and rric, bnt is independent of p and p^. Similarly, 
we denote by still a generic positive constant to address the dependence on p. In addition, 
a prodnct space (W)” of vector fnnctions is still denoted by X, for example, a vector fnnction 
u E {HP is denoted hj u E H^ with norm \\u\\h^ := (X]fc=i ■ 

2.1. Incompressible case 

Onr first main resnlt is concerned with the instability/stability for the linearized incompress¬ 
ible magnetic RT problem fll.4p - fll.6p and reads as follows. 

Theorem 2.1. Denote the critical number by 



m, 


r .= snp 

Dhi 


g f p'wjdx 
XofjdiWpx 


(i = 1 or 3). 


( 2 . 1 ) 


Assume that the density profile p satisfies (II.2p . Then, uiq is a threshold of M for instability 
and stability of the problem fll.4p - fll.6p in the following sense: 

(1) If \m\ < m^, then the equilibrium state {p,0,M) of the problem fll.4p - fll.6p is unstable. 
That is, there is an unstable solution 


{g, u, N, q) := e^fi-p'us/^, u, mdiU, q) 

to fll.4p - fll.6p . where {u,q) E fl .4,^) X H^ solves the boundary value problem: 

j X^pu = ApAu —'\/{Aq + XQm'^diUi) + Xom’^dfu + gp'usCs, divu = 0, 

\ u\dn = 0 

with a finite growth rate A > 0 satisfying 

= snp Srj{w,K) = £„{u, A), 


( 2 . 2 ) 


(2.3) 
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where 


Moreover, 


S„{w,A) := J gp'w^dx — Xom^ J \diWfdx — Ap J |Vwfdx. 

b(^)||L2, \\{ui,U 2 ){t)\\L 2 , \\u 3 {t)\\L 2 , 

||(A^1, A^2)(^)||l2 and ||A^3(t)||i2 co as t ^ oo. 


(2.4) 


(2) If \m\ > rriQ and the initial data qq^uq^Nq satisfy 

Mo G n H^, {go, No, diNo) G and divA^o = 0, 

then there is a unique global solution {g,u,N) G x C^{'R'^,H'^) x to 

the problem fll.4l) - fll.6p satisfying the following stability estimates: for any t > 0, 

\\{g,u,N)t{t)\\l 2 + \\{u,diu){t)\\l 2 +p [ WurllHidr <C\\{go,diUo,iaAuo,diNo)\\l 2 , (2.5) 


Ht)\\h + 


L2 


{u,diu){T)dT 

\\{g,N){t)\\l2<C\\{go,uo,No,Qo)\\l2, 
/^l|VM(t)||i2 < C\\{go, diUo, pAuo, No, diNo)\\l2, 


+ A \\u{T)\\HidT <C\\{uo,Qo)\\l2, 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


where Qo := XomdiNo — gog^s- Moreover, there exists a couple (^?oo,^oo) G If, such that 

j{XomNoo ■ dip + ggooq:> 3 )dx = 0 for any p = {pi, p 2 , q:> 3 ) ^ H^, (2.9) 

||2(^) - 2 oo||l 2, ||M(t)||j:^i, \\ut{t)\\L 2 , ||iV(t) - iVoolUs0 ast^oo. (2.10) 

Remark 2.1. It should be remarked that we have omitted to write the restriction w 7^ 0 in 
fl2.ll) in order to make the denominator sense. We mention the dehnition fl2.1l) makes sense, and 
niQ G (0, +00). In fact, recalling the condition fll.2l) for p, there exists a function w G Hf, such 
that 

J p'wldx > 0, 

which guarantees the validity of fl2.ip and implies tuq > 0. On the other hand, we have 

WwiWi^ < cn\\djWi\\L 2 for w e Hq, (2.11) 

which implies < +00. Here cq is a positive constant depending on the diameter of O only. 

Remark 2.2. For the linearized nonhomogeneous incompressible viscous fluid problem, it is 
shown in 20|, Theorem 1.1] that a RT equilibrium state is always unstable. However, Theorem 
[Q shows that in the presence of a sufficiently large magnetic field M, the linear stability can 
be maintained due to the stabilizing effect of the magnetic field. In addition, when the domain 
is symmetric about the plane xi = X 3 (e.g., a ball centered at the original point) and p' is a 
constant, then 


i gp 

mQ = — sup 

-^0 w&H}, 


J w^dx 


(* = 1,3), 


/ \diw\‘^dx 

and obviously, mQ = mf,. This shows that the horizontal magnetic field has the same stabilizing 
effect as the vertical one, a new phenomenon found in this article. 
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Remark 2.3. Theorem 12.11 still holds for a horizontally periodic domain with hnite height (i.e., 
n := (27rLT)^ x (—/,/), I G (0,+oo)). In this case, we can assert that 

ml = Me, (2.12) 

which will be shown in Section [5] directly by using the dehnitions (12.Ih and (II.7p . In addition, 
we can see from fl2.ip that ml = +cxd for the horizontally periodic domain. This means that 
the horizonal magnetic held has no stabilizing effect (due to the horizontally periodic boundary 
condition). 

Remark 2.4. We briehy describe the idea in the derivation of (12.4p . It follows from (12.2p 
that M G Mq n satishes = £„{u,A) > 0; moreover, u G C°’‘'(r2) for some constant ^ G 
(0, 1) by Sobolev’s embedding theorem. Thus, one has p'us 7 ^ 0, hs 7 ^ 0 and diUs 7 ^ 0, and 
consequently, ||p(t)||i2, ||M 3 (t)||L 2 , ||^ 3 (^)|| L2 —)• 00 as t —>■ 00. Recalling that divw = 0 and 
u\dn = 0, we get 7^ 0 ^ind \diU\\‘^ + \diU 2 \‘^ 7^ 0 immediately, which yields \\{ui,U 2 )(t)\\L 2 

and ||(W, A^2)(^)||l2 —)■ 00 as t — )■ cx). Such process is consistent with the phenomenon of the 
magnetic RT instability: gravity hrst drives the 3-th component of velocity unstable. Then, the 
instability of the 3-th component of the velocity further results in instability of the density and 
horizontal velocity. Finally, the instability of the velocity leads to instability of the magnetic 
held through the induction equation. 


Remark 2.5. To our best knowledge, it is still open whether there exists a small smooth solution 
of an initial boundary value problem for the equations of a homogeneous incompressible viscous 
MHD how of zero resistivity without external forces in a bounded domain (i.e., p is a constant 
and p = 0 in fll.ip b Therefore, at present we can not establish a nonlinear stability result for the 
more complicated problem fll.3l) - fll.6p . On the other hand, for the nonlinear instability problem, 
even if we have a local-in-time existence result for fll.3l) - fll.6p . there still are some difficulties 
induced by boundary in establis hing the nonlinear instability from the linear instaffility by a 

we only 


induced by boundary m establis hing 

classical bootstrap argument BQ. Due to the same reason, in a previous work 
considered a horizontally periodic domain with inhnite height. We mention that for the Cauchy 
problem, one does have certain existence results on global small smooth solutions 
large solutions [^. 


30 and local 


Remark 2.6. In this remark we give some extended results. 

• If D is a C^+^-smooth bounded domain and p G C'*’’''^(r2)jthen {u,q) G x by 

the classical regularity theory on the Stokes problem (see [9|, Theorem IV.6.1]). 

• If m = m^, then we only have the following stability estimate of Ut, which can be observed 
from fl2.18p . 

\Wt\\\2+pl \\ur\\]jidT <C\\{Qo,diUo,p/Auo,diNo)\\\2. (2.13) 

Jo 

• If p is further in C^(r2), then we can deduce that dj{p'u^) ^ 0 for j = 1,2,3, and hence 
||^jT(^)||l2 —)■ C )0 as f —)■ 00. Moreover, based on 02.61) . we can deduce from 01.6p that the 
solution in the linear stability result satishes ||9jp(t)||^2 < C\\{go, djgo,uo, Qo)\\l 2 , whence, 
lim^^oo \\djgit) — c?jPoo||L2 = 0, also see the derivation of 03.19p . 

• Obviously, Theorem 12.11 still holds for i = 2 and M = (0,m,0). 





























Next, we sketch the main idea of the proof of Theorem 12.11 to explain how to utilize the 
critical number, and the detailed proof will be given in Section 121 We start with the proof idea 
of the linear instability. As in 2^, we make the following ansatz of growing mode solutions to 
the linearized problem: 


g{x,t) = g{x)e^^, u{x,t) = u{x)e^^, = g(x)e^‘, N{x,t) = N{x)e^^ for some A > 0. 


..Ai 


At 




Substituting this ansatz into fll.bp . we get 
( Kg + p'us = 0, 

< Apu + V(g + XomNi) = pAv + XomdiN — gge^, divu = 0, (2-14) 

[ AN = mdiU, diviV = 0, 


and then eliminating p and N by using the hrst and third equations, we arrive at the time- 
invariant boundary problem fl2.2p for u and q. Then we apply a modihed variational method as 
in to construct a solutions of the eigenvalue problem fl2.2p . This idea was used probably hrst 
by Guo and Tice to deal with an ODE problem arising in the construction of unstable linear 
solutions 


16j . In view of the basic idea of the modihed variational method, we modify 


as 


follows. 


r apu = spAu -|- gp'u^e^ + Vp -|- Xom'^dfu, divu = 0, 

I hlan = 0, 


(2.15) 


where s > 0 is a parameter, p := —{sq -|- Xorn?diUi) and a := a(s) depends on s. The system 
fl2.15p satishes the identity: a{s)J{u) = £„{u,s), where 


£cr{u, s) := E^{u) — sp / iVupda: and E^{u) := / gg'n^dx — Xom? / \diu\^(lx. 


Thus, using the variational method and the classical regularity theory on the Stokes problem, we 
hnd that u G O Aa, and a satishes fl2.15l) by maximizing 


q;(s) = sup £„{w,s) G M. 

WGAa 


( 2 . 16 ) 


Moreover, in view of the dehnition of q:(s), we can infer that a(s) G oo) is nonincreasing. 

On the other hand, recalling the instability condition \m\ < ttiq, we conclude that E„{wo) > 0 for 
some Wq G and thus there exists a hnite interval (0, ©) on which a{s) > 0 and lim^^e <a(s) = 
0. Employing a hxed-point argument, we immediately see that there is a A satisfying 

A = \/a(A) = / sup £^{w,A) G (0, 6). (2.17) 

V weAa 


Consequently, we can construct a nontrivial solution u G to fl2.2p with a positive A dehned 
by fl2.17p . and therefore the linear instability follows. 

For the linear stability, the key observation is that any solution to the linearized magnetic 
RT problem fll.4p - fll.6p satishes 


j p\ut\‘^dx 


E^{u(t)) 2p 



\VUr\‘^dxdT = Jo, 


(2.18) 
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where Jq := / (pl^iP + — gp'u^) |t=oclx. On the other hand, recalling the dehnition of 

the critical number, we get 



\diu\^dx^ 


whence, 

- E^{u) > Ao[m^ - [rn^'^f] j \diu\^<lx. 

Hence fl2.5l) follows from fl2.19p . fll.hp i . fll.61) -:. and the condition \m\ > m\.. 
Furthermore, similarly to the derivation of fl2.18l] . we hnd that 



dx + 2p / |VM(r)|^dxdr 


p\u\ -E^[ M(r)dr 


= /p|Mo|^da: + 2 Qoudxdr < / p|Mo|^dx + 2 ||Qo||l2 



M(r)dr 


L2 


(2.19) 


( 2 . 20 ) 


which, together with Cauchy-Schwarz’s inequality, gives fl2.6p immediately. Finally, with these 
stability estimates in hand, we can deduce fl2.7p - fl2.10p by standard energy estimates and asymp¬ 
totic analysis. 


2.2. Compressible case 

Based the ideas in the proof of Theorem 12.11 we can further establish the instability/stability 
of the linearized Parker problem fll.14p - fll.16p by a more careful analysis for compressibility. 
Unfortunately, we can not give a critical number of the magnetic held for fll.14p - fll.16p as for 
the incompressible problem fll.4p - fll.6p . since the horizontal steady magnetic held is vertically 
stratihed. Now, we introduce another version of instability/stability criterion for fll.14p - fll.16p . 

First, substituting the following ansatz of growing mode solutions into fll.161) . 

q{x) = Q{x)e^^, u{x) = u{x)e^^, N{x) = N{x)e^^, 


we get 

' Ag + div(p-u) = 0, 

Apu + V{p'{p)g + AomciVi) = pAv + /ioVdivh -h AoiVgM' Xom^diN - ggcs, 

AN = rricdiu — u^M'^ — Mcdivu, 

, diviV = 0, 

and then eliminating p and by using fl2.21l) i and 02.211) 9. we arrive at a time-independent 
boundary value problem problem for u = {ui,U 2 ,U 3 )'. 

{ A?pu = gp'usCs + V[p'(p)div(p-u)] -|- gpdivues -|- Xomcirricdlu — Mcddvdiu) 

+V[XQmc{jncd 2 U 2 + mcd^u^ + m'^u^)] + ApAu + ApoVdiv-u, (2.22) 

h|ao = 0. 
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Multiplying fl2.22p by u, integrating the resulting equations and using the condition fll.lip . 
we infer that 


J pluf dx = J [gp'ul — p'(p)p|divMf + Xomc{mcdiu ■ u — rricUidwdiu) 

— Xomc{mcd2U2 + mcd^u^ + m'ji3)d\vu 

+ gpuzdivu — p\p)p'uzdivu]dx ~ ^ J + /Uo|div-u|^)dx 
= J {gp'ul — p' {p)p\divu\‘^ — Xoml[\diu\‘^ + {d2U2 + 93M3)divu — ^lUidivu] 
+ {gp — Xorricm'^ — p'(p)p')u3divu}dx ~ ^ J + /^o|div-up)dx 
= j [gp'ul - p'(p)p|divM|^ - Xoml{\diU2\^ + \diu^\^ + \diUi - divu]^) 

+ 2gpu^ddvu]dx — A / {p\Vu\^ + /io|divM|^)da;. 


Dehning 


Eciu) := j gp'uldx + j 2gpd\vuu2,dx — J p'(p)p\divu\‘^ 
+X(iml{\diU2\‘^ + IdiUsl'^ + \d2U2 + dx, 

V{u) =; / {p\Vu\‘^ + po\divu\‘^)dx, 


we have the relation: 

AV(m) = S^{u, A) := E^{u) - AV{u). (2.23) 

Recalling the proof of linear instability in Theorem 12.11 the instability condition \m\ < thq 
guarantees that Erj{w) > 0 for some w G H{. This gives the existence of a positive growth-rate 
A by applying the modihed variational method. Similarly, in view of (I2.23p . we shall impose the 
following condition: 

there exists a function w G H^, such that Edw) > 0, (2.24) 

which is actually a necessary condition for the existence of a positive growth-rate A. Of course, 
we shall show that fl2.24p indeed is the instability condition for the linearized Parker problem in 
Subsection 14.11 

Now we turn to a sufficiently condition for the linear stability. In the proof of the linear 
stability in Theorem 12.11 the key step is to deduce the relation fl2.19p . from which and the 
condition |mc| > tuq we obtain the desired stability estimates immediately. Dehne 

^ _ Edw) _ 

I Ao(| 5 iW 2 P + + 152^2 + 

where we have omitted to write the restriction that w G should make the denominator and 
the square root operation sense. Hence, if Cr < 0, we hnd that 

Edu) < AoCr J (|9iM2p + + |^2M2 + dsUsl^dx for any u G ifg. (2.25) 
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Then we can deduce from fl2.25p the stability estimate on the velocity of solutions to the linearized 
Parker problem, the detailed derivation will be presented in Subsection 14.21 In addition, it is 
easy to check that Cr > 0 is equivalent to fl2.24|) : moreover Cr 7 ^ + 00 . 

Now, we sum up the previous discussions and state the instability and stability results for 
the linearized Parker problem fll.l4|) - fll.l6p . 

Theorem 2.2. Assume that p and rric satisfy fll.2p and fll.lOp . and let Me = (?7ic,0,0). 

(1) //Cr > 0, then the equilibrium state (p, 0, Me) of the problem fll.l4l) - fll.l6p is unstable, that 
is, there exists an unstable solution 

(p, u, N) := e^*(—div(p-u)/A, u, {mcdiu — u^M'^ — Medivu)/A) 

of fll.14p - fll.16l) . where u G //^ D A solves the boundary value problem fl 2 . 22 p with a finite 
growth rate A > 0 satisfying 

= sup A) = £c{u,A). (2.26) 

w&A 


Moreover, 

||(Mi,M 2 )(t)||L 2 , \\u 3 {t)\\L 2 , ||(A^ 1 , A^ 2 )(t)||L 2 and \\N 3 {t)\\L 2 ^ 00 as t ^ 00 . 

In addition, ||p(t)||L 2 -A 00 as t —)■ 00 provided p' > 0. 

(2) If Cr < 0, and the initial data go,uo,No satisfy 

Uo e n H\ go eH\ {No, diNo, VN^) e L'^ x L'^ x L'^ and divNo = 0, 

then there exists a unique global solution {g, u, N) e C°(R+, H^) x C'°(R+, //^) x C°(R+, 
of fll.i 4 p - fll.l 6 l) satisfying the following stability estimates: 


\\{ 0 ,u,N)t\\l 2 + IIdiva)(t) 11^2 + p / ||'U,-||^idT < CIo, 


(2.27) 


Ht)\\h + 


{u, diu, divM)(r)dr 


+ p ||'a(r)||^idr < C'dKrto, Po)|Il2 +/q), ( 2 . 28 ) 


L2 


||(p, A^)(t )||^2 < C{\\{go,Uo, Nq, Po)\\‘l 2 + /q), 
p||Vw(t)||i2 <C(||(po,wo,iVo,Po)||i2+/o), 


(2.29) 

(2.30) 


where Pq = V{p'{p)go + XoitVcN^) — XoAf^M^ — XomcdiNo + gogoo, and Iq is defined in fl4.14l) 
and can be bounded from above by 


C[bo||Hi + (p + Po)||Am||l2 + \\{diu°,d2U2,d3ul,Nl,diNo,VN^)\\L2]. 
Furthermore, there is a couple {goo,Noo) € L^, such that 

j [XottIcNoc ■ dif) + ggoofJz - XoN^mA^i - (p'(p)poo + AomcA^“)div^/>]da: 

for any ip = (^/’i,^> 2 , V’s) e Hq, 

ll^'(^) - ^’ooIIlz, ||M(t)||j^i, ||Mt||L 2 and ||//(f) -//oo||l 2 -)■ 0 as t-)■ 00 , 
where (j = 1,3^ denotes the j-th component of N°°. 


(2.31) 


(2.32) 

(2.33) 
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Remark 2.7. We can choose a sufficiently large Imd, such that Cj. < 0. In fact, we dehne 


I + “^gpuisdivui — p'{p)p\diYw\^)dx 

wgh^ V I Ao(|9lW2p + + \d 2 W 2 + 93W3|2)dx' 

By Remark l2.ll we see that n G (0,+ 00 ). In view of fll.lOp . for given p, p, g and Aq, we can 
choose rric satisfying infa;gn |^c| > h. Then, it is easy to verify that such rric satishes Cr < 0. 

Remark 2.8. We give some extended results similar to those in Remark 12.61 

• If is a C^^^-smooth bounded domain and p G then u G by the classical 

regularity theory on elliptic equations (see jlOl. Theorem 4.11]). 

• Any linear solution of fll.14p - fll.16p satishes the following estimate 

\\VpUt\\l2 + 2 p [ ||VMr||i2dr < Jo, 

Jo 

provided Edw) < 0 for any w E Hi (i.e., Cr < 0). 


3. Proof of Theorem 12.11 

3.1. Linear instability 

In this Subsection we give the detailed proof of the linear instability in Theorem 12.11 To 
begin with, we show that a maximizer of fl2.16p exists and that the corresponding Euler-Lagrange 
equations are equivalent to fl2.15p . 

Proposition 3.1. Assume that the density profile p satisfies the first two conditions in fll.2p . 
then for any but fixed s > 0, the following assertions are valid. 

(1) £a{w,s) achieves its supremum on Aa-. 

(2) Let u be a maximizer and a satisfy fl2.16p . then there exists a corresponding pressure field 
p associated with u, such that the couple (u, p) satisfies the boundary value problem fl2.15p . 
Moreover, {u,p) G x . 


Proof. (1) Let {wnyidi C Aa be a maximizing sequence, then £a{wn,s) is bounded. Recalling 
that 

£a{wn,s) = j[gp'iw^fi - Xom‘^\diWn\‘^ - sp\Vwfi‘^]dx, 

we can easily see that Wn is bounded in Hfi Here wf denotes the third component of Wn- So, 
there exists a w E Aa and a subsequence (still denoted by Wn for simplicity), such that Wn ^ w 
weakly in Hi and strongly in L^. Moreover, by the lower semi-continuity, one has 


sup £a{w,s) = limsup £’o-(tan, s) 

wGAc n—>oo 


= lim 

n^oo 


p'{w'^fdx— liminf / {Xo7rr\diWn\ + sp\Vwn\)dx 


<£a{w,s) < sup £a{w,s), 

WGAct 


which shows that £a{w, s) achieves its supremum on Aa- 
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(2) To show the second assertion, we notice that since S^{w, s) and J{w) are homogeneous of 
degree 2, fl2.16p is equivalent to 

£Aw,s) 

a = sup ■ (3.1) 

w&Hl J 

For any r G M and ip G H^, we take v{t) := u + Tip, where h is a maximizer. Then (13.Ih implies 

s) - J(h(r)) < 0. 

If we set I{t) = S„{v{t),s) — A‘^J{v{t)), then we see that /(r) G C'^(M), /(r) < 0 for all r G M 
and /(O) = 0. This implies 

dJ(r) 


dr 


= 0 . 


T = 0 


Hence, a direct computation leads to 


a pu ■ ipdx = — sp Vu : Vipdx — Xom^ / diU ■ diipdx + / gp'u^e^ ■ ipdx, 


which implies that u G is a weak solution to fl2.15p . Thus, by virtue of scaling the variable 
Xi, it follows from the classical regularity theory on the Stokes problem that there are constants 
Cl dependent of the domain fl, s, p, Aq and m, and C 2 dependent of ci, g, a and p, such that 

I|w||h2 + ||Vp||l 2 < cilK^fp'hsCs - apu)\\L2 < C2, 

where p G FfMs the corresponding pressure held associated to u. This completes the proof. □ 

Next, we want to show that there is a hxed point s = A > 0 such that ^/a{K) = A. To this 
end, we hrst give some properties of a(s) as a function of s > 0 . 

Proposition 3.2. Under the assumptions of Proposition lS^ the function a(s) defined on (0, oo) 
enjoys the following properties: 

(1) Q!(s) G C|^’j!(0, cxd) is nonincreasing. 

(2) If \m\ < niQ, then there are constants c^, C 4 > 0 which depend on g, p, p, Aq and m, such 
that 

a(s) > C 3 — SC 4 . (3.2) 

Proof. (1) Let {wn}'^=i C Aa be a maximizing sequence of sup^g_ 4 ^ Sj) for j = I and 2. 
Then 


tt(si) > limsup -^i) — liminf - 52 ) = Q;(s 2 ) for any 0 < Si < S 2 < cxo. 

n—>-oo n^oo 

Hence a(s) is nonincreasing on (0, 00 ). Next we use this fact to show the continuity of a(s). 

Let I := [a, b] C (0, 00 ) be a bounded interval. In view of the monotonicity of a(s), we know 
that 

|a(s)| < max{|a(a)|,p IIp'/pII^oo} < C )0 for any s G/. (3.3) 

On the other hand, for any s G /, there exists a maximizing sequence C Aa of 

sup^g_ 4 ^ £a-(uj, s), such that 

|a(s)-£:^( m;^,s)| < 1. (3.4) 
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Making use of fl3.3p and 03.41) . we infer that 


0 </x J [Vwpdx + J\diw\‘^dx 

=» [ pWax - 

S J s 

^ l + niax{|a(a)|,^||pVp||^^} ^ ^ 
~ a a 


:= K, 


where denotes the third component of w^. Thus, for Sj E I {j = 1, 2 ), we further hnd that 
q;(si) = hmsup£^o-(tt'^\ Si) < limsup'^ 2 ) + hki ~'S 2 I hmsup / |Vta^^|^dx 

n^oo n^oo n—>-oo J 

<Q!(s2) + K\si — S2|. 


(3.5) 


Reversing the role of the indices 1 and 2 in the derivation of the inequality 03.5p . we obtain the 
same boundedness with the indices switched. Therefore, we deduce that 

|a(si) - Q!(s2)| < K\si - S2I, 

which yields a{s) G (^100(0, cxd). 

(2) We turn to the proof 03.2p . Noting that m < rriQ, we can deduce from the dehnition of 
tjiq that there is a n G H^, such that 


gp'v^dx — Xom^ / \divpdx > 0 . 


Thus, one has 


/ \ ^ / \ ^crid^i s) 

a(s) = sup t^[w,s) = sup ———— 

weAa •U1&HI J \w) 


> 


8a{v, s) f gp'v^dx — Xovn? f \div\^dx p f | Vnpdx 


J(v) 


fpv^dx 


— s- 


f pv'^dx 


C 3 — SC 4 


for two positive constants C 3 := C 3 {g,p,m) and C 4 := C 4 (/i,p). This completes the proof of 
Proposition 13.21 □ 


Next we show that there exists a pair of functions {u, q) satisfying 02.21) with a growth rate 
A > 0. Let 

G := sup{s I a(r) > 0 for any r G (0, s)}. 

By virtue of Proposition 13.21 G > 0; moreover, q;(s) > 0 for any s < G (in addition, we can 
further show that a(s) strictly decreases on (0,G)). Since a{s) = sup^g_ 4 ^ To-(tc, s) < 00 , using 
the monotonicity of a(s), we see that 


lim q;(s) exists and the limit is a positve constant. (3.6) 


On the other hand, by virtue of Poincare’s inequality, there is a constant C 5 dependent of g, p 
and O, such that 


9 


p'w^dx < C 5 


|Vtc|^da; for any w E A. 
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Thus, if s > CsZ/i, then 


9 


p'wldx — sp 


iVwpdx < 0 for any w & A, 


which implies that 


a(s) < 0 for any s > c^/p. 


Hence © < oo. Moreover, 


lim a(s) = 0. 

s^6 


(3.7) 


Now, exploiting (I3.6p . (13.7p and the continuity of a{s) on (0,©), we hnd by a hxed-point 
argument on (0,©) that there is a unique A G (0,©) satisfying 02.171) . Thus, by virtue of 
Proposition 13.11 there is a solution [u, q) G x to the problem 02.2p with A constructed in 
02.17p . where q := —(p + Xom^diUi)/A. We conclude the following proposition, which, together 
with Remark 12.41 yields the linear instability in Theorem 12.11 


Proposition 3.3. Assume that the density profile p satisfies 01.2p and m < m^. Then there 
exists a pair of functions {it, q) G (77^ fl Aa) x 77^ which solves the boundary value problem 
02.21) with a finite growth rate A > 0 satisfying 02.31) . In particular, g := —p'u^/A G 77^ and 
N := mdiii/A G Moreover, {g,it,N,q) satisfies 02.14p . 


Remark 3.1. It is easy to see that Proposition 13.31 still holds for p being a function of three 
variables {xi,X 2 ,x^). However, if 

Vp = -pge^, 

then p has to be a function of the single variable x^. This is the reason why we only consider 
that p is a single variable of x^ in Theorem 12.11 


3.2. Linear stability 

Before proving the linear stability in Theorem 12.11 we shall establish the local well-posedness 
of the linearized magnetic RT problem 01.4p - 01.6p . which can be shown by an iterative method. 
Next, we briefly describe how to show it for the reader’s convenience. 

Let T* G (0,1), It* ■= (0, T*), K > 0 and 

UK = {ue C\lT*,H‘^)nC\lT*,Hl) I < 77}, 

where T* and K will be hxed later. Given v G Uk, we consider the following linear problem: 

r put + Vp = pAu + f, divM = 0, 

1 u\dQ = 0, u\t=o = Uq, 

where 

n 3 (r)dr + po, N = m f div{T)dT + Nq. (3.9) 

Jo 

Obviously, divA^ = 0, since divAg = 0. Using fl3.9p . we have {f,ft) G C^{It*,LJ), and thus the 
problem fl3.8l) possesses a unique solution u G L°°{It* , H"^) with a unique associated pressure 


/ = XomdiN - gges, 9 = - p' 

Jo 
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p G LP') satisfying Jpdx = 0. Moreover, u and p enjoy the following estimates (referring 


to 


j^, Lemma 5] or 22, Lemma 5.2]) 


Il“i|lc0(/*,L2) + + ll“illL2(/j,.,/fi) + \\P\\\^{I,J.* ,m) 

< C'm[II“o|Ih2 + II/o|Il 2 + (1 + ^*)||/|lLoo(/^,^i2) + 7'*||/i|||oo(/^,_2.2)] (3.10) 

< ^(ll^»o||i 2 + ||mo ||^2 + \\diNQ \\\2 + K‘^T*) for any t e It*, 


where the constant C only depends on /i, p, Aq, g, m and the domain f2. Using the classical 
regularity theory on the Stokes problem, we have 


\\u{t2) -u{ti)\\H2 + \\p{t2) -p{ti)\\Hi < C^{\\Ut{t2) -Mt(tl)||L2 + ||/(t2) - /(tl)||L2), 

which implies that u G C^{It*, H"^) and p G C^{It*, H^). Now, taking K > (5(||po|li2 + ||wo ||^2 + 
||9iiVo||^2 + 1) and T* < 1/K'^, we arrive at m G Uk- 

Considering the above results, we can construct a function sequence satisfying 


r pdtUn+i + Vpn+i = pAw„+i + XomdiNn - Qngez, divM„+i = 0, 

1 ^n+l\dfl O 5 ^n+l|t=0 ^0 


and 





\\9tUn\\cO(I^^,L^) + IK^nl 

2 - + 

WdtUnWh^T 

where 

pt 




Qn = - p’Vs 

Jo 

(r)dr + po, 

Nn = m 

Let 

^n+1 ^n+1 Qn Qn 

~ 0n—l, Hn ~ 

1 

II 


we have 


pdtUn+l + Vpn+1 = pAun+1 + XomdiN^ - Pnge^, diVM„+i = 0, 
^n+ll^O 0, Un-\-l\t=0 0 


In view of the estimate (13.101) . we get 

\\dtUn+l\\ ,L'2) + 

+ + ||Pn+l||L°°(/y.,rfi) < CT*\\Un\\Loo(^j^^^H2y 


Hence, we see that {un}’^=i, {dtUn}’^=i and {pn}'^=i are a Cauchy sequence in 

L\lT*,H^)nL^{lT*,L^) and 


for sufficiently small T, respectively. Thus we obtain the limit functions m, p. It is easy to verify 
that the limit functions u, p are a unique solution to the following problem 

{ pUt + Vp = p/S.u + XomdiN — gge^, divw = 0, 

P = - Jo p'u 3 (r)dT + po, A = m f* diu(T)dT + Nq, (3-11) 

u\dn = 0, u\t=o = uo- 


Furthermore, 


ueC°{lT*,H^), uteC%lT*,L^), 



0 and p e H^)- 


17 





Obviously, {g, u, N) constructed above also uniquely solves the linearized problem fll.4p - fll.6p with 
an associated pressure q = p—X^mNi. Moreover, (p, iV, diN) G due to (po, Nq, diNo) G 

To get a global solution in Theorem 12.11 it suffices to deduce the global estimates for ||p(t)|| 2 , 2 , 
||M(t)||ji ^2 and \\diN{t)\\i 2 . To begin with, we derive the energy equality fl2.18p . In view of the 
regularity of {g,u,N), we can deduce from fl3.11|) i that for a.e. f > 0, 


lA 

2ffi 


p\ut\ dx =< putt, ut >= / {Xom^dju • Ut + gp'U 2 ,dtU 2 )dx -p \Vut\ dx, 


where < •, • > denotes the dual product between the spaces and and represents 
the dual space of H^, please refer to [^, Remark 6]. On the other hand. 


lA 

2ffi 


{Xom^\diu\‘^ - gp'ul) = J ■ diUt - gp'usdtUs) dx 

= - {gp'usdtus + Xom^diU ■ ut) dx. 


Putting the previous two equalities together, we conclude 

^ J + Xom^\diu\‘^ - gp'ul) dx + 2/i J [Vutl^dx = 0, 

which yields fl^TTSP . 

Then, using the inequality fl2.19p . we further infer from fl2.18p that 

J {plutp + Ao[m^ — (my^] dx + 2/i J J |VnT-pdxdr < Jq. 

Recalling fl2.1ip and Poincare’s inequality, we obtain 


WutWh + \\{u,diu)\\l2 +p lln^ll^idr < CJo- 


(3.12) 


Applying 03.121) to the hrst and third equations in 01.6p . we hnd that 


\\{g,u,N)t\\l 2 + p\\{u,diu)\\l 2 + p / WurWH^dr < CJo- 


(3.13) 


On the other hand, we have 


Jo < \\{0o,diUo, pAuo,diNo)\\l^2. 


(3.14) 


In fact, using the second condition in 01.21) and Cauchy-Schwarz’s inequality, we obtain from 
Ol.bD o that 


p\ur\^{T)dx = / {pAu + XomdiN — gge^) ■ Ur{T)dx 


inf p 


<C{\\giT)\\i2 + \\pAuir)\\i2 + II Wr)||i.) + ^||n.(r 


Il2, 


which implies 


|Mr(r)||i 2 < C'(||p(r)||i 2 + \\pAu{t )\\12 + \\diN{T)\\l 2 ) for any r > 0. 
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Let r —)■ 0, then 


||Mt|t=o||i2 < (711 (^0) ^i-^o)|lL2) 


and one gets fl3.14p . Putting fl3.13p and 03.141) together, we obtain the stability estimate 02.Sp . 

Noting that diN = m d‘fu{T)(iT + 9jiVo, thus, we use the classical regularity theory on the 
Stokes problem to infer that 

\\u{t)\\H^ + ||g(t)||r/i <Cf,{\\ut{t)\\L 2 + IIAom^iA^ - gge^Wi^) 

<C^^(^\\{go,diUo,Auo,diNo)\\L2 + j^ ||M(r)||^^2dr^ , 
which, together with Grownwall’s inequality, yields 


\\u{t)\\m < C^\\{go,diUo,Auo,diNo)\\L2 (l + . 

With the global estimate 03.15P and 


(3.15) 


IKe,aiiv)(i)|| <C IKeo,9iiVo)llL> + / ll«(T)||H=dT 


in hand, we immediately get the global solution {g,u,N) with an associated pressure g by a 
continuity argument based on the local well-posedness result. Moreover, the global solution 
satishes the stability estimate 02.51) . 

Now, we proceed to deriving the estimates 02.6p - 02.8l) . Firstly, Ol.bp o yields 

put + V(g + XomNi) = pAu + Xorn^dt / M(r)dr + p' M 3 (r)drge 3 + Qq. 

Jo Jo 

Consequently, O2.20p follows. Noting that u{T)dT\g^ = 0, similarly to the derivation of 03.12p . 
we can obtain 02.61) . Hence, exploiting 01.6|) i . we have 


||2WIIl2 <||2o||l 2 + / 2rdr 

Jo 

<I|2oIIl2 + IIp'IIl 

Similarly, using Ol.bp i. one obtains 

||A^(t)||L2 <||fVo||L2+ / NAt 


L2 

t 


M3(r)dr 


(3.16) 


< C\\{go,Uo,Qo)\\L2. 


L2 


L2 


<||A^o||l2 + \m\ 


9j-u(r)dr 


(3.17) 


^ C\\{uo, Nq, Qo)\\l^. 


L2 


Therefore, the estimate 02.7p follows from 03.16P and 03.17p . Using 02.5|) - 02.7p . we deduce from 
(fL6l) 9 that 

p\\Vu\\l 2 + \\VpUt\\l 2 = - {XomN ■ diU + ggu:i)dx 


<fi'||2||L2||M||L2 + XQm\\N\\L2\\diU\\L2 
^C'll {goi diUOi pAUf), Nq, 5 jiVo)||L 2 ) 


(3.18) 
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which yields fl2.8p . Hence, to complete the proof of the linear stability in Theorem 12.11 it remains 
to show the asymptotic stability of {g,u,N) in fl2.1Up . 

By j^, Theorem 1.68], Sobolev’s inequality, the estimates on ||m|| in fl2.5p and fl2.7p . 

we infer that 


— ||(m, VM)||i 2 dr < 2 / {\\u\\l 2 \\Ut\\h-^ + \\'^u\\L^\\VUr\\H-^)dT 


<C \\{u,Vu,Ur,VUr)\\l2dT 


<C \\{u,Ur)\\HidT < C\\{go,diUo, fxAuo,diNo)\\l2. 


Hence, ||■u(^)||i^l G H^^’^(0, cxd), which, together with (I3.18p . implies ||M(f)||Hi and ||Mi(t)||2,2 —)■ 0 
as t —)■ cxD. Using fl3.16p and (13.171) . we see that there are two measurable functions Poo and 
and a time sequence {tn}'^=i C M"*", such that 


(p, N){tn) —)■ (poo 5 Nao) weakly in L‘^ as n ^ oo. 


On the other hand, we have that for any t G M"*", 


|| 2 (^) - 2 oo||l 2 < liminf ||p(f) - p(f„)||i 2 < liminf / UprlUadr 


(3.19) 


<IIp1 


WusWi^dr < C'||mo||l2 


and 


||A^(t) - A^oo||l 2 < / ||A^r||L 2 dr < m / ||9jM(r)||L2dr < 0 ||mo||l 2 . 


Thus, II (p — Poo, N — iVoo)(t)||L 2 —)■ 0 as f —>■ cx). Finally, multiplying (11.6p by u, we obtain 


— / pUt ■ (fdx = / (/iVw : V</9 + XomN ■ dip + ppe 3 )da;. 


which, together with (12.101) . implies (12.91) . The proof of the linear stability in Theorem 12.11 is 
complete. 


4. Proof of Theorem 12.21 

In this section, we adapt the basic ideas in the proof of Theorem 12.11 to show Theorem 12.21 
Due to the compressibility, the proof of Theorem 12.21 will be more complicated than that of 
Theorem 12.11 

4.I. Linear instability 

We still apply a modified variational method to construct a solutions of the boundary value 
problem (I2.22p . so we modify (12.221) as follows. 

{ apu = spAu + s/ioVdiv-u + gpdivues + gp'u^es + V[p'(p)div(p-u)] 

+Aomc(mcC?i'U — Mcdivcii-u) + 'V[Xomc{mcd2U2 + + m'-us)], (4.1) 

n|(9o = 0, 
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where a := a(s) depends on s. Then, the standard energy functional for fl4.ip is given by 


Tc(m, s) = Ec{u) - si/( m) 

with an associated admissible set A. Thus, we hnd an a by maximizing 


a := sup s). 
w^A 


(4.2) 


Obviously, sup^g_ 4 Tc(tc, s) < oo for any s > 0. Next we show the existence of a maximizer for 

(O. 


Proposition 4.1. Assume that the density profile p satisfies the first two conditions in fll.2p . 
then for any but fixed s > 0, the following assertions are valid. 

(1) £c{w,s) achieves its supremum on A. 

(2) Let u be a maximizer and a defined by fl4.2l) . then u G is a weak solution to the boundary 
value problem gHi, 

(3) lfa>0, then the maximizer u further satisfies 


its 7^ 0, mcdius ^ 0, (4.3) 

Im'hg + mc{d2U2 + + |mc5iM2p 7^ 0, (4.4) 

uj + Ll^O, (4.5) 

where Ut denotes the i-th component of it. In addition, 

div(p-u) 7 ^ 0, provided p' > 0. (4.6) 


Proof. (1) Let <Z A he a. maximizing sequence. Using Cauchy-Schwarz’s inequality, 

one sees that 


£c{Wn,s) =EfiWn) - sV{Wn) 


< 

< 


gp'\wf\‘^dx + J ( 25 fpdivt(;„t (;3 — p'(p)p\divwn\‘^) dx — sp 
gp'\wf\‘^dx + J ^ ^ dx — sp J \'Vwn\'^dx, 


\VWn\‘^dx 


which yields 


sp 


\VWn\‘^dx + SfiWn, s) < 


gp'\wf I'^dx + 


9^p\w3 

p'{p) 


2 

-dx. 


We easily see from the above inequality that Edwn, s) is bounded, and consequently, Wn is 
bounded in Hq. So, there exists a w & A and a subsequence (still denoted by Wn for simplicity), 
such that Wn ^ w weakly in Hq and strongly in L^. Moreover, by the lower semi-continuity, one 
has 


sup £c{w, s) = limsupTc(wn5 s) 

w£A n—>-00 


= lim / igp'\w'f\^dx + 2gpdi'vWnW'^) dx— 

n—>-oo J ^ n—^oo 

+\oml{\diW 2 \‘^ + \diw^\'^ + \d 2 W 2 + dsw'^l'^) + s{p\Vwn\‘^ + /io|divM;„|^)] dx 

<£c{'dj,s) < sup s), 

w£A 


j [p\p)p\divWn 
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which shows that £ciw, s) achieves its supremum on A. 

(2) To show the second assertion, we write 04.21) as follows 


a = sup 

w&Hl 


£c{w,s) 

J{w) 


For any r G M and '0 G Hq, we take v{t) \= u + T'lp. Then, 04.Th implies 

Ec{v{t),s) - \^J{v{t)) < 0. 

Let /(r) = Ec{v{t), s) — A^J(i)(r)), then 


(4.7) 


d/(r) 


dr 


T = 0 


0 . 


Hence, a direct computation leads to 

a J pu ■ 'ipdx = j[{gp'u'ie'i + ^fpdivues) • ip + gpdi'vipe2, ■ u\dx 

- j {\Qml[diU2diip2 + diu^idiip^ + {d2U2 + d3U3){d2ip2 + 

+ p'{p)pdivudivip}dx ~ ^ J (aVm ; Vip + podivMdiv' 0 )da;. 

Using the condition Ol.llh . we change the above weak form as follows 

a J pu ■ ipdx = — J[{spo + p'{p)p)divu + \oml{d2U2 + d3U3)]divipdx 
— sp J Vm : V'^da; + Xoml J{divudiipi — diu ■ diijj)dx 

+ / [{gp'uses + gpdivue3 + V{p'{p)p'u3) + XqV{ rricm'^us)] • ipdx, 


(4.8) 


which implies that F is a weak solution to (14.11) . 

(3) Next, we turn to the proof of (14. 3 D - (14. 6 j) by contradiction. 

By the second assertion, we know the maximizer u E A satisfies 04.Sp . thus a = Sc{.u,s). 
Suppose U 3 = 0 , then a = £c{u, s) < 0 due to J p\u\‘^dx = 1 , which contradict with the condition 
a > 0. Hence M 3 7 ^ 0. 

Suppose |m(.M 3 + mc{d 2 U 2 + + |mc 5 iM 2 p = 0, then 


m'^us + mc{d2U2 + 83^3) = 0 


(4.9) 


and mcdiU 2 = 0. Since rric > 0, we get diU 2 = 0. Recalling that U 2 \dn = 0, we obtain U 2 = 0 
by a Lagrangian formula [33|, Section 1.3.5.1], which, tougher with 04.91) . implies d 3 {mcU 3 ) = 0. 
We immediately see M 3 = 0, which contradicts with M 3 7 ^ 0. Hence 04.4p holds. Similarly, we can 
show mdiU 3 7 ^ 0 . 
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Suppose that u\ + U 2 = Q oi div(pM) = 0, then 


0 <a= J ^gp'ul + 2 gpu 3 d 3 U 3 - (p'(p)p + Aom^) - Aom^l^iMal^jda; 

-s j (pIVusI^ + dx 

= - j [{P'{p)p + ^oml)\d3U3\^ + \oml\diU3\^]dx 
-s (pIVmsP+ po|53M3n dx < 0, 


or 


0 < a = J ^gp'ul + {2gpu3 — p'(p)pdiYu)divu 

- \oml{\diU 2 \‘^ + IdiUsl"^ + \d 2 U 2 + dsu^l'^^dx - s J (p| Vup + poldiv-up) dx 

= - j [ 9 P'ul+p'ip)p\divu\‘^ + Aom^d^iMsP + + \d 2 U 2 + ^gUs^jda; 

— s y (p|Vu|^ + po|div-ud) dx < 0, 

which is a contradiction. Therefore, fl4.5p and fl4.6l) hold. This completes the proof. 


□ 


Proposition 4.2. Let u G he a weak solution of the boundary value problem fl4.ip . then 
u G H"^. 

Proof. Firstly, we write fl4.ip as follows. 

{sp + spo + p'{p)p)dlui + spdlui + spdlui 
+(spo + p'{p)p)did2U2 + (s/io + p'{p)p)did3U3 = apui + gpdiu^, 

{sp + Xoml)dfu2 + {sp + spo + p'{p)p + Aom^)9|-U2 + spdlu2 
+ {spQ + p'{p)p)did2Ui + (s/io + Xoml + p'{p)p)8283113 = apu2 + 

{sp + Xoml)8fu3 + SP82U3 + 83[{spo + p'{p)p)8iUi 
+{spo + XquiI + p’{p)p)82U2 + {sp + spQ + Xoml + p'{p)p)83U3] 

= apU3 - gd\v{pu) + 83{gpU3), 

P\dn = 0, 

Let / = (opui + gp 8 iU 3 , apU 2 + pp^aUs, - 5 'div(pM) + 83 {gpU 3 )Y and be the 

matrix of coefficients of the linear elliptic equations 04.101) . then O4.10p can be written as 

-8a{Aff8yUj) = fi, 

where we have used the Einstein convention of summing over repeated indices, and the non-zero 
coefficients are 

322 — — 4^^ — 4^^ — ‘ill 

^11 ~ ^11 ~ ^22 ~ ^33 ~ 

AI2 = = All = Alsi = spo + p'(p)p, 

^22 = ^33 = sp + XomI, All = All = spo + p'(p)p + Aom^, 

= sp + spo p'(p)p. All = sp + spo + p{p)p Xoml, 

A33 = sp + spo p'{p)p + Xoml. 
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Noting that, for any f] E 

+ (s/^ + P'(p)p)(6hi + 6h2 + ^sVsY 

+ Aom2(6h2 + + Kml^Kv2 + vl) > 

hence satishes the strong elliptic condition. On the other hand, dfl is of class C"^, G 
0°’^(r2) and / G Thus, if we apply jl^. Theorem 4.11] to the weak form fl4.8p . we get u G if^. 
□ 


Similarly to Proposition 13.21 we have the following properties of a(s) as a function of s > 0. 


Proposition 4.3. Under the assumptions of Proposition^^ the function a(s) defined on (0, oo) 
enjoys the following properties: 


(1) Q!(s) G cxd) is nonincreasing. 

(2) //Cr > 0, then there are constants Cq, Cy > 0 which depend on g, p, p, po, Aq, p and rUc, 
such that 

a{s) > cq — scj. 


Proof. The monotonicity and the second assert obviously hold by directly following the proof 
of Proposition [221 while the absolute continuity can be established by modifying the proof of the 
hrst assertion in Proposition 13.21 Here we give the proof of absolute continuity for the reader’s 
convenience. 

Let / := [a,b] C (0, cxd) be a bounded interval. For any w G Al, by Cauchy-Schwarz’s 
inequality, 

Sc{w,s) < / {gp'wl + 2gpdivww3)dx — / p'{p)p\diYw\^dx 


<9 


+ 


9P 


p'{p) 


Hence, from the monotonicity of Q!(s) we get 


|q!(s)| < max < \a{a)\, g 


+ 


9P 


p'{p) 


:= L < oo for any s G /. 


(4,11) 


On the other hand, for any s G /, there exists a maximizing sequence {tc*} C of sup^g^ £c(w, s), 
such that 

|a(s)-£’c«,s)| < 1. (4.12) 

Making use of fl4.11|) and fl4.12p . we infer from the dehnition of Sc{w, s) that 

0<y (/i|Vtcp +/ioldivwp) dx 

<- [ [c/p'Ksr + -p'(p)MvOdiv<]da; “ 


/1 9 

1 

p' 

+ 

L°° 

9P 


1 + 2L 

'' 1 
a a 

p 

p'{p) 

L°°- 

“ a 


Thus, for Sj G / (i = 1, 2), we have 

a(si) < a{s 2 ) + (1 + 2L)|si - S 2 |/a, 
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and further infer that 


|a(si) - a(s 2 )| < (1 + 2L)\si - S 2 \/a, 


which yields a(s) G <^ 100 ( 0 , C)o). This completes the proof of Proposition 14.31 


□ 


With Proposition 14.31 in hand, we can directly follow the proof of fl2.17p to deduce that there 
is a unique A > 0 satisfying 


A = \/ a(A) = ^ I sup Sc{w, A) > 0. 

wGA 


(4.13) 


Thus, by virtue of Propositions 14.11 and 14.21 there is a solution u G H'^ to the boundary value 
problem fl2.22p with A constructed in 04.131) . where u satishes 04.311 - 04.61) . Thus, we conclude the 
following proposition, which gives the linear instability in Theorem 12.21 

Proposition 4.4. Assume that the density profile p satisfies 01.21) and Cr > 0. Then there 
exists a u E D A which solves the boundary value problem 02 . 22 |) with a finite growth 
rate A > 0 satisfying 02.26p . Moreover, u satisfies 04.311 - 04.61) . In particular, let {p,N) : = 
(—div(p'u),mc 5 i'u —usM' —Mcdiv-u)/A, then {p,u,N) G L'^ x H'^ x L'^ solves 02.211) . iVf+ iV| 7^ 0 
and A^3 7^ 0. In addition, p 7^ 0 provided p' > 0. 

f. 2 . Linear stability 

Firstly, following the process of the iterative method as in Subsection 13.21 we can show that 
there exists a unique local-in-time solution (p, u, N) of the linearized Parker problem 01.141) - 
01.16p . where we should use the classical regularity theory on elliptic equations instead of that 
for the Stokes problem (see 04.191) later). Moreover, the solution satishes the following regularity 

p G C\It*,H^), (iV,aiiV, ViVi) G C\It*,L^), 
u G C\It*,H‘^), ut G C\lT-^,L^)nL^ilT-^,H^), 

where T* G (0,1) depends on (||po||^i + ||mo ||^2 + IK-^S) ^ 1 -^ 0 , VA^{’)||^ 2 )- Next, we deduce some 
global-in-time estimates on ||p(t)||_H-i, ||M(t)||/f2 and ||(Ai"3, c^iA^, VA^i)(t)||^2- 

In view of the regularity of (p, u, N), we deduce from 03.1ip i that for a.e. f > 0, 

PUtt,Ut > 

= j[{p'{p)Qt + Aomc9Wi)divMt -F Xom'^dtNsdtUi 

+ Xom^diNt ■ Ut - QtgdtUsldx - j{p\Vut\‘^ + po\dwut\'^)dx := Li + L 2 , 

where < •, • > denotes the dual product between the spaces H~^ and Hq, and H~^ denotes the 
dual space of Hq. On the other hand, using 01.16p i . 01.16p o. 01.lip and partial integrations, we 
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have 


Li = J [—p'{p)div{pu) + XorricijncdiUi — m'paz — mcdwu)\dwutdx 

+ j[\om^^{d‘lu ■ Ut — didivudtUi) + gdiY{pu)dtU'^dx 

= j {gp'usdtus - (p'(p)p' + \omcm'^)usdivut + gpdivudtUs - p'{p)pdivudwut 

— Xoml[diu ■ diUt + (divw — diUi)divut — divudtdiUi]dx 

= j {gp'u^dtUs + gpdtidivuus) - p'{p)pdwudivut 

- Xoml[diU 2 dtdiU 2 + diu^dtdiUs + {d 2 U 2 + d 3 U 3 ){d 2 U 2 + d 3 U 3 )t\dx 

=-^EJu(t)). 

2dt ^ ^ ” 

Combining the above two equalities, we get 


1 d 
2dt 


which yields 


p\ut\^dx - Ec{u{t)) + / {p\Vut\^ + po\divut\‘^)dx = 0, 


{p\Vur\‘^ + /io|divM^p)da:dr 


p\ut\‘^dx - Ec{u{t)) + 2 



(4.14) 



= J (pl^tTL^o “ •= ^ 0 - 

Keeping in mind that ||VdivMo||L 2 < ||Amo||l 2 (see [HI, Corollary 9.10]) and (I2.1ip . we easily see 
that Jo can be bounded from above by (I2.3ip . 

Recalling fl2.25p . we have 

[ p\ut\‘^dx - AoCr J{\diU 2 \^ + + \d 2 U 2 + d 3 U 3 \^)dx 

{p\Vur\^ + /io|divMT-|^)dxdr < Jq. 

Since Cr < 0, we obtain by exploiting Poincare’s inequality and fl2.1ip that 

ll“t(^)llL2 + \\{u2,U3,diU2,diU3,d2U2 + d3U3){t)\\\2 

+ / (/i||Mr||^i + ho||div-u.r||i 2 )dr < CJo, for any t e 
Jo 

On the other hand, fl4.14p can be rewritten as follows 

[ [plwtT +p'(p)p|divu|^ + Xoml{\diU2\'^ + + \d 2 U 2 + ^sUap)] dx 

(/xIVutP + /io|divM.rndxdr = Jq + / ( 25 fpdiv-uM 3 + gp'uf) dx. 


(4.15) 



Thus, using Cauchy-Schwarz’s inequality and fl4.15p . we obtain 

lldivull^s < C/o, 


(4.16) 
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which, together with fl4.15p and fl2.1ip . yields 


\\{ui,dim)\\l 2 <CIo (4.17) 

Putting fl4.15l) - fl4.17l) . fll.lhp i and fll.lhp o together, we get the stability estimate fl2.27p . 

Now, we turn to the estimation of ||M(t)||ji/2. To this end, we use fll.lhp i and fll.lGp -i to rewrite 
fll.lbp o as follows. 

fiAu + /xoVdivM = put + V {p'{p)p + XorricNi) — XqN^M'^ — XorricdiN + gge^, (4.18) 


where 


Q=— div(pu){T)dT + Qo and N = / {rricdiu — u^M'^ — Mcdivu){T)dT + Nq. 

Jo Jo 

The viscosity term in fl4.18p defines a strong elliptic operator on u, thus 
\\u{t)\\H 2 <\\pUt + V (p'(p)p + Xom^Ni) - AoiVgM^ - Xom^diN + 






||pobi + ||(iV30,aiiVo,ViV{',M,)llL2+ / \HT)\\H^dT 


(4.19) 


||pobi + ||(iV3^5i^o,ViV°)|U2 + Jo+ / ||M(r)||^,2dr 


which, together with Grownwall’s inequality, gives 

Wumn^ < + \\{Nl,d,No,VN^„)\\L 2 + Io) (l + . 


(4.20) 


Here denotes a generic positive constant depending on G, p, pq and the other known 

physical parameters. With the help of fl4.20p and 

I|£>(*)IIh. + IKA's.SiA'. V]Vi)(()|li> sc (ikolli,. + IKAfstSiiVo, V]V»)|U. + f ||«(T)||„=dT) , 


we immediately obtain a global solution {g,u,N) by a continuity argument based on the local 
well-posedness result. Moreover, the global solution satisfies the stability estimate fl2.27p . 

We proceed to deriving the estimates fl2.28p - fl2.30p . Firstly, we get from fl4.18l) that 


p'(p)div ( p / M(r)dr ) + X^nic / {nicdiUi — m[us — mcdivM)(r)dr 


put + V 


= pAu + poVdivu + Xorricdi / {rricdiu — u^M'^ — McdivM)(r)dr 

Jo 

+ Xorricdi / M3(r)drM' + ^fdiv (p M(r)dr ) 63 + Pq- 


Similarly to the identity fl2.18p . we have 


p\u\‘^dx — Ec i / M(r)dr ) + 2 



{p\Vu\‘^ + /io|divMp)(r)da;dr 


(4.21) 


= / p\uo\‘^dx + 2 Po M(r)drda;. 
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Noting that /g K(r)dr|g^ = 0, we obtain (I2.28p from (14.211) by following the derivation of (I2.27p . 
Utilizing (ll.lbp i . we hnd that 


||2(^)IIl2 <I|2o||l2 + 
= I|2o||l2 + 


div(pM)(r)dr 


div ( p / M(r)dr 


L2 


(4.22) 


< C'dI (po, Wo, .Po)||l 2 + Jo). 


L2 


Similarly, using fll.lbp q. one obtains 


l|iVWIlL2 <||a^oIIl2 + 


{rricdiu — u^M'^ — McdivM)(r)dr 


L2 


(4.23) 


<C{\\{No,uo,Po)\\l^ + Io)- 


Combining (Km with we arrive at (Km . 

Now, making use of fl2.28l) . fl2.29l) and Cauchy-Schwarz’s inequality, we infer from fll.lbll o that 

p\\Vu\\l2 + po\\divu\\l2 

= — J[put ■ u — {p'{p)q + AomcA^i)divM — Xqui'^N^Ui + XorricN ■ diu + ggusldx 

< C{\\{g,N3,Ut)\\LAH\L^ + ||(p, A^i)||L2||divM||2.2 + \\N\\L2\\diu\\L2) 

< C(||( 20 ,Wo,iVo,Po)||L 2 + /o), 

which implies (12.3011 . In addition, following the proof of (12.911 and (12.1011 . and using the stability 
estimates, we obtain the asymptotic behaviors (12.3211 and (I2.33p . The proof of linear stability 
results in Theorem 12.21 is complete. 


5. Additional results 


5.1. Critical number of horizontally periodic domains 

In this subsection we prove the equality (12.1211 in Remark 12.31 Obviously, it suffices to show 
the following conclusion. 

Proposition 5.1. Let L > 0, Z > 0, O := (27rLT)^ x {—1,1), p := p'{x^) G L 

f p'w^(x)dx 


a = sup 


w{x)&Hi / \d3w{x)\‘^dx 


and b = 


sup I , 

Pix3)eH^i-l,l) J_i\i’'{x3)\‘^dX3 


f_l p'\^/J{x3)\^dX3 


then a = b. 


Proof. Let W 3 {f,X 3 ) be the horizontal Fourier transform of W 3 {x) G Hf, i.e.. 


W3(,^,X3) 


'(27rLT)2 


W3{x',X3)e ^^''^dx', 


where x' = {xi,X 2 ) and f = (^i,6), then dsWs = dsths. We denote 'ip{f,X 3 ) := 'ipi{f,X 3 ) + 
i'ip 2 {^,X 3 ) := W 3 {f,X 3 ), where and V ’2 are real functions. By the Fubini and Parseval theorems 
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(see jl2|, Proposition 3.1.16]), we have 

1 


^ / (p'l{p/>o}+p'l{p'<o})|w3(C,a:3)rda;3 


«G(L- 


1 

'4^2 


X] / P'(^i(^,a;3) + V’2(C,a;3))da:3 


(5.1) 


i&{L- 


and 


1 


A 

1-1 


\d3W3{x)\^dx = I + \d3^2{^,X3)\‘^)dx3, (5.2) 


?e(L-iz)2 

where l{p>o} and l{p<o} denote the characteristic functions. Noting that -01 £ by the 

dehnition of b, we have 

b[ \d3^IJi{^,X3)\‘^dX3> [ p'\^jJi{^,X3)\‘^dX3, 


'-I 


'-I 


Thus, using fl5.ip - fl5.2p . we immediately deduce that 


^ ^ Eg6(L-iZ)2 I-l P'\'iPi^,X3)\‘^dX3 _ fp'\w3{x)\‘^dx ^ / p'\w3{x)\^dx 


E.g(L-izp I-l \d 3 H^,X 3 Wdx 3 I | 53 n^ 3 (a:)Pda; / |a 3 n;(a;)| 2 da; ’ 


Next we turn to the proof of a > 6. We choose a maximizing sequence C ifg(— /,Z) 


Hence a < b 
Next we 
of b to see at 


I-iP'\i’j{x3)?dx3 
I-l\'^j{x3WdX3 


b = lim 


(5.3) 


For each xjjj, we can construct an approximate sequence {'0m}m=i ^ 0 satisfying 

in /) as m —)■ oo. Now, we denote 


w = (O,'0(„(a;3) cos(?7.L ^X 2 ),nL Vm(a^ 3 ) sin(?7,L ^^ 2 )), 


then w{x) G and 

J p'wl{x)dx 3 


f p'w3(x)dx 


f |93M;(a;)|2da; f(jd3W2(x)j^ + |53t(;3(a;)|2)dx 
_ (nL"i )2 Y P'^mi^s) sin^(nL-^X2)dxidx2dx3 

cos(nL-ia;2)P + |nL-Vm(a^3) sm{nL-'^X2W)dxidx2dx3 
I-ip''<Plix3)dx3 


L‘^n-^ /!; \ij'^{x3)\^dx3 + Y \ij'^{x3)\‘^dx3 


for sufficiently large m, 


whence. 


. ,. ,. I-iP'i^iix3)dx3 f_ip'i^j{x3)dx3 

a > hm hm --. 

m^ocn^oo ^ 2^-2 \^'^{x3)\^dX3 + | ^’^( 0 : 3 ) | 2 d(r 3 \ij'j{x3)\^dX3 


which, together with (15.3p . yields a > b. The proof is complete. 


□ 
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5.2. Sharp growth rate of solutions to the linearized problems 

In this section we show that A dehned by fl2.3p resp. fl2.26p is the sharp growth rate for any 
solution of the linearized problem fll.4p - fll.6p resp. fll.14p - fll.16p . We shall exploit the energy 
estimates as in 1^, [lOl, l20| to show that is indeed the sharp growth rate for {g,u,N) in 
X X L^-norm. 

Proposition 5.2. (i) Let {g,u,N) solve the linearized magnetic RT problem fll.4p - fll.6p with an 
associated pressure q. Then for any t >0, 

||2(t)||x<C^/.e"^‘(ll2o||x + ll(Vwo,Awo,aWo)||i 2 ), X = or H\ 

\\u{t)\\jji + \\ut{t)\\l 2 + [ ||VM(r)||^ 2 dr < Vmq, Amq, 9Wo)|li2, 

Jo 

||iV(f)|U2 < Vko, Ano,iVo,aWo)||i2, 

where the constant only depends on p and A. 

(a) Let (p, u, N) solve the linearized Parker problem fll.14p - fll.16l) . Then for any t >0, 

MtWm + htimh + f l|V«(r)||i 2 dr < + hi 

Jo 

11 ( 2 , A)(f)||i 2 < Ao)||i 2 + holl^i +/o), 

where the constant only depends on /i, po A. 


(5.4) 

(5.5) 


Proof. We prove only the second assertion and the hrst assertion can be shown in the same 
manner. Let {g,u,N) be a solution of fll.l4l) - fll.l6l) . then {g,u,N) satishes the identity fl4.14p . 
In view of fl2.26l) . we have 

Ec{u) < hfj{u) + A / (/ilVhp + /io|div-up)da;, 


which, combined with fl4.14l) . results in 


J p\ut\‘^dx + 2 J J (p\Xut-\‘^ + po\divUr\‘^)dxdT 
< h + J p\u\‘^dx + A J {p\Vu\‘^ + po\divu\‘^)dx. 

Using Newton-Leibniz’s formula and Cauchy-Schwarz’s inequality, we hnd that 
A{p\\Vu{t)\\l 2 + po\\divu{t)\\l 2 ) 


— Kn + 2A 



0 JQ 


T ^x^^^-Uj^x^^^-UjdxdT + podivM,-divM j dxdr 

1<*J<3 

pt 


<K„ 


(A'IIVUtIIIs +;iol|divMT||ii)dT +/ (ai||Vm(t)|||i +/ici||divu(r)|||j)dT, 


(5.6) 


(5.7) 


where Kq := A(p||V mo ||^2 + ho||cliv«o||^ 2 )- Thus, we infer by fl5.6p - fl5.7p that 
1 


A 


\\y/put{t)\\L2 + R\\Xu{t)\\L2 + po\\dwu{t)\\L2 


< A||\/pM(t)||i2+ 2A / (/i||VM(r)||i2+po||divM(r)||i2)dr + 


h + ‘JKq 
A 


(5.8) 
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Recalling that 


A^||\/pM(t)|li 2 = 2A j pu{t)-ut{t)dx<\\y/put{t)\\l 2 +A‘^\\^u{t)\\l 2 , 
we further deduce from fIS.Sp the differential inequality: 


■^\\Vpu{t)\\l2 + p\\Vu{t)\\l2 + po\\diyu{t)\\l2 


< 2A 


\\Vpu{t)\\l2+ / (/i||VM(r)||^2+/io||divM(r)||i2)dr 


+ 


In + 2Kn 


A 


Applying Gronwall’s inequality jssi, Lemma 1.2] to the above inequality, one concludes 


\\Vpu{t)\\l2+ / (/x||VM(r)||22+;io||divn(r)||i2)dr < 
Jo 

which, together with fIS.Sp . yields 


WVp'^oWl^ + 


2 , (-^0 + 2A’o) 


2A2 


,2 At 


(5.9) 


j\\VpUt{t)\\h + hl|VM(t)||i2 + /xo||divn(t)||2 2 


< 2 


MWf^oWL^ + 


2 I (-^0 + 2iLo) 


2A 


2At I -^0 + 2iLo 
^ A • 


Thus fl5.4p follows from the two estimates above. Finally, using fll.lbp i . fll.lbp o and fl5.9p . we 
hnd that ^ 

II(2,^)(^)||l2 <||(2o,^o)||l2 + [ ||(p, A^)r||L2dr 


<II(2o,^o)||l 2 + C" / ||M(r)||Hidr 
Jo 

No)\\l^ + ||Mo||rti + V%), 

and get fl5.5p . The proof is complete. 


□ 
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